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Recently, we have predicted [Phys. Rev. Lett. 97, 226102 (2006)] that a single ethylene molecule 
can form stable complexes with light transition metals (TM) such as Ti and the resulting TM„- 
ethylene complex can absorb up to ~12 and 14 wt % hydrogen for n=l and 2, respectively. Here we 
extend this study to include a large number of other metals and different isomeric structures. We 
obtained interesting results for light metals such as Li. The ethylene molecule is able to complex 
with two Li atoms with a binding energy of 0.7 eV/Li which then binds up to two H2 molecules per 
Li with a binding energy of 0.24 eV/H2 and absorption capacity of 16 wt %, a record high value 
reported so far. The stability of the proposed metal-ethylene complexes was tested by extensive 
calculations such as normal-mode analysis, finite temperature first-principles molecular dynamics 
(MD) simulations, and reaction path calculations. The phonon and MD simulations indicate that 
the proposed structures are stable up to 500 K. The reaction path calculations indicate about 
1 eV activation barrier for the TM2-ethylene complex to transform into a possible lower energy 
configuration where the ethylene molecule is dissociated. Importantly, no matter which isometric 
configuration the TM2-ethylene complex possesses, the TM atoms are able to bind multiple hydrogen 
molecules with suitable binding energy for room temperature storage. These results suggest that 
co-deposition of ethylene with a suitable precursor of TM or Li into nanopores of light-weight host 
materials may be a very promising route to discovering new materials with high-capacity hydrogen 
absorption properties. 



PACS numbers: 68.43.Bc, 81.07.-b,84.60.Ve 

I. INTRODUCTION 

The success of future hydrogen and fuel-cell technolo- 
gies is critically dependent upon the discovery of new 
materials that can store a large amount of hydrogen 
at ambient conditions Recently, from quantum me- 
chanical calculations we found that the C=C bond in 
a single ethylene molecule, similar to and carbon 
nanotubes^^^^, can form a stable complex with tran- 
sition metals (TM) such as Ti. 9 . The resulting TM 2 - 
ethylene complex attracts up to ten hydrogen molecules 
via the Dewar-Kubas interaction 1 ^ reaching a gravimet- 
ric storage capacity of ^14 weight-percent (wt The 
interaction between hydrogen molecules and transition 
metals lies between chemi- and physi- sorption, with a 
binding energy of ~0.4 eV/H2 compatible with room 
temperature desorption/absorption. 

Different from metal decorated Cgo/nanotubes, metal- 
C2H4 complexes are actually existing structures and have 
been actively studied in the past several decades, with 
the major goal being to understand the catalytic mech- 
anisms and processes of metals. Experimental spectro- 
scopic data on various complexes, such as Li, Mg, Al 
and TMs complexed with C2H4, widely exist in the 
literatureiii 1 ^^ 3 -^. These complexes were typically syn- 
thesized by direct reaction of metal atoms with C2H4/Ar 
in the gas phase. Early theoretical studie a 14 ' 15 ' 16 ' 17 ' 18 
showed that the metal-C2H4 binding mechanisms could 
be either electrostatic (e.g., C2H4-AI), or Dewar-Chatt- 
Duncanson bonding (e.g., most C2H4-TMS). Despite the 



extensive studies on the metal-C2H4 complexes, their 
ability to absorb H2 was not realized and investigated 
until our recent work 9 -. 

Here we extend our earlier work 9 - and present a de- 
tailed theoretical study of the hydrogen absorption on a 
large number of metal-C2H4 complexes, including TMs 
and the alkali metal Li. We organize the paper as fol- 
lows. In the next section, we describe the computa- 
tional methodology. In Sec. Ill, we discuss C2H4M com- 
plexes, various isomers of C2H4M2 complexes and present 
the metal binding energies, zero-temperature lattice dy- 
namics of these complexes and their hydrogen absorp- 
tion properties (including the H2 binding energies and 
maximum number of H2 that the complex can absorb). 
In Sec. IV, we discuss the possible reaction paths (i.e., 
minimum energy paths) and the activation energies (i.e. 
barriers) between various isomers of C2H4Ti2 complexes. 
We also discuss an interesting catalytic effect of Ti, sim- 
ilar to the "spillover effect" , where a molecularly bound 
H2 molecule is first dissociated over Ti and then one of 
the H atoms is bonded to carbon, forming a CH3 group. 
The resulting molecule is isostructural to an "ethanol" 
molecule and thus called "titanol" . The titanol molecule 
is also able to absorb up to five H2 as molecules with a 
binding energy of ^0.4 eV/H2 and provide another in- 
teresting possibility for high-capacity hydrogen storage 
materials. In Sec.V, we present high-temperature first- 
principles molecular dynamics (MD) studies on selected 
structures. Due to the small system size, we are able to 
carry out MD simulations up to 10 ps. We show that the 



proposed complex structures are quite stable and exhibit 
constructive desorption upon heating without destroy- 
ing the underlying complex. Our concluding remarks are 
presented in Sec. VI. 



II. DETAILS OF CALCULATIONS 

Our first-principles energy calculations were done 
within density functional theory using Vanderbilt- 
type ultra-soft pseudopotentials with Perdew-Zunger 
exchange correlation, as implemented in the PWscf 
packagei 1 - Single molecular complexes have been treated 
in a supercell of 20x20x20 A with T k-point and a cut- 
off energy of 408 eV. The structures are optimized until 
the maximum force allowed on each atom is less than 
0.01 eV/ A for both spin-paired and spin-relaxed cases. 
The reaction path calculations were carried out using the 
Nudged Elastic Band (NEB) methoc" We used a to- 
tal of 21 images between the reactant and the product, 
which were fully optimized during the NEB calculations. 
The MD simulations were carried out within the micro- 
canonical ensemble (NVE) starting with the optimized 
structure and random initial atom velocities'*^ More 
details of the MD calculations are given in Sec.V. 



III. STRUCTURAL, ELECTRONIC, AND 
DYNAMICAL PROPERTIES OF C 2 H 4 M n AND 
C 2 H 4 M„-H :c COMPLEXES 

We start by examining various possible configurations 
of C2H 4 M rl complexes and their corresponding H2 ab- 
sorption properties. We consider both transition metals 
and light metal Li, and focus on n=l and n=2 cases. 
Complexes with n >2 are less attractive for hydrogen 
storage due to potentially lower capacities and thus are 
not discussed here and should be avoided in the synthe- 
ses. 

When one metal atom binds to the ethylene molecule, 
the configuration shown in Fig. 1(a) is the most ener- 
getically favorable one, where the metal atom forms a 
symmetric bridge "bond" with the C=C bond of ethy- 
lene. When two metal atoms bind to C2H4, the complex 
may adopt several possible configurations. In our ini- 
tial study£, we focused on the sandwich structure (Fig. 
1(b)). Here we consider two additional isomeric struc- 
tures: dimer-par (Fig. 1(c)) and dimer-perp (Fig. 1(e)). 
In the sandwich configuration, each M atom is closer to 
one of the carbon atoms, leading to two different M-C 
"bonds". Note that for most transition metals (e.g., Ti), 
there is no classical chemical covalent bonding between 
the metal atom and carbon atom. The calculated bond- 
population is found to be nearly zero for these metals. 
The slight shift of the metal atoms towards different C 
atoms only results in a minute contribution of the M-C 
covalent-like bond to the overall binding. In just a few 
cases (e.g., Fe), the metal and carbon atom are bonded 
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(a) C 2 H 4 -M (b) C 2 H 4 -M 2 (sandwich) 




FIG. 1: (Color online) Various configurations of C2H4M n 
(n—1 and 2) complexes considered in this study, (a) C2H4 
complexed with one metal atom, (b)-(e) C2H4 complexed 
with two metal atoms with different metal binding sites. Note 
that the bond-stick model is only used for clarity and should 
not be considered as an implication of the chemical covalent 
bonding between those atoms. For most metals, there is no 
classical chemical covalent bonding between the metal and 
carbon atoms. For a few metals (e.g., Fe), the complexes pos- 
sess a structure, where M and C are bonded more traditionally 
by covalent bonding, as shown in (d). (f) C2H4M2 complex 
with dissociated C=C bond. Large, medium and small balls 
represent M, C and H atoms, respectively. 



more traditionally by a covalent bond, as shown in Fig. 
1(d). For this reason, we generally specify these C2H4M„ 
structures as "complexes" instead of " molecules" . 

The binding mechanism of the C2H4TM complex has 
been discussed in detail in our previous work£. Essen- 
tially, the bonding orbital for the TM- atoms and C2H4 
results from the hybridization of the lowest-unoccupied 
molecular orbital (LUMO) of the ethylene molecule and 
the TM-tf orbitals, in accord with Dewar coordination. 
For Li, the binding mechanism is, however, different. In 
Fig. 2, we show the electronic density of states of C2H4 
molecule, Li atom and C2H4Li complex. Projection anal- 
ysis of the states indicates that the electron in the 2s 
state of Li is divided into two halves that are transferred 
to the LUMO of C2H4 and the 2p of the Li atom, respec- 
tively. Then the 2p orbital of Li and the LUMO of C 2 H 4 
are hybridized for the binding of Li on the C2H4. From 
the isosurfaces of the molecular orbitals (also shown in 
Fig. 2), it is clear that the molecular orbital of C2H4Li 
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TABLE I: The metal-C2H4 binding energies (in eV/M-atom) with respect to atomic and bulk energies of various metals, and 
the average H2 binding energies (in eV/Hb) on C2H4M for various absorption configurations (see Fig. 3). The maximum 
number of H2 molecules bonded to each metal is also shown. 



Property /M 


Li 


Sc 


Ti 


V 


Cr 


Mn 


Fe 


Co 


Ni 


Cu Zn 


Zr 


Mo 


W 


Pd 


Pt 


E B (M-atomic) 


0.32 


1.39 


1.45 


0.94 


0.18 


0.51 


0.92 


1.39 


0.91 


0.80 none 


1.91 


1.02 


1.71 


1.95 


2.52 


E B (M-bulk) 


-1.41 


-2.72 


-3.68 


-4.30 


-3.44 


-3.06 


-1.64 


-2.43 


-1.99 


-2.86 


-4.23 


-5.19 


-6.65 


-1.86 


-2.82 


E B (per H 2 ), MH 2 




0.96 


1.16 


1.00 


0.01 


0.59 


1.01 


0.94 


1.13 


0.19 


1.90 


0.85 


1.80 


0.83 


1.33 


E B (per H 2 ), M+H 2 


0.29 


0.02 


0.31 


0.46 


0.45 












0.35 


0.49 


0.59 


0.64 




max H 2 /M 


2 


5 


5 


5 


5 


5 


5 


3 


2 


2 


5 


5 


5 


2 


2 


E B (per H 2 ), M+2H 2 


0.28 
















0.42 


0.33 








0.27 


0.25 


E B (per H 2 ), MH 2 +3H 2 




0.40 


0.54 


0.66 


0.34 


0.24 


0.31 








0.78 


0.61 


0.85 






E B (per H 2 ), M+5H 2 




0.28 


0.46 


0.53 


0.21 


0.18 


0.34 








0.54 


0.64 


0.79 







TABLE II: The metal-C2H4 binding energies (in eV/M-atom) of three isomeric C2H4M2 configurations (see Fig. 1), with 
respect to atomic and dimer energies of various metals. 



Property /M 


Li 


Sc 


Ti 


V 


Cr 


Mn 


Fe 


Co 


Ni 


Cu 


Zn 


Zr 


Mo 


W Pd 


Pt 


Eb (M-atomic), Sandwich 


0.69 


1.39 


1.47 


1.21 


0.05 


0.37 


0.83 


1.30 


0.70 


1.41 


none 


1.69 


0.37 


1.18 1.56 


1.78 


Eb (M-atomic), Dimer-par 


0.54 


1.77 


2.02 


1.62 


0.10 


0.64 


1.65 


1.63 


1.09 


1.34 


none 


2.66 


2.20 


3.26 1.88 


2.61 


Eb (M-atomic), Dimer-perp 


0.61 


1.72 


2.12 


1.97 


0.21 


0.51 


1.22 


1.50 


1.00 


1.24 


none 


2.70 


2.10 


2.41 1.39 


1.53 


Eb (M-dimer), Sandwich 


0.20 


0.58 


0.17 


-0.21 


0.79 


0.33 


-0.36 


-3.45 


-0.01 


0.17 




-0.10 


-1.71 


-1.30 0.75 


0.20 


Eb (M-dimer), Dimer-par 


0.05 


0.96 


0.72 


0.20 


0.84 


0.60 


0.46 


-3.12 


0.38 


0.10 




0.87 


0.12 


0.78 1.08 


1.03 


Eb (M-dimer), Dimer-perp 


0.12 


0.91 


0.82 


0.54 


0.95 


0.47 


0.03 


-3.25 


0.29 


0.00 




0.91 


0.02 


-0.07 0.58 


-0.05 



near the zero-energy (i.e., the Fermi Energy), is a su- 
perposition of the LUMO of C2H4 and the p-orbital of 
the Li atom. Also note that the occupied orbital of the 
C2H 4 Li complex at around -4eV is about the same as that 
the highest occupied molecular orbital (HOMO) of bare 
C2H4, except that there is a hole in the upper portion of 
the orbital due to the Li-ion. The bond analysis does not 
show any covalent bonding between C and Li atoms. For 
C2H4Li2, we observed also a binding mechanism similar 
to that of C2H4LL 

The metal binding energies on ethylene are summa- 
rized in Table I and Table II for one metal and two metal 
complexes, respectively. They are calculated by subtract- 
ing the equilibrium total energy Et of the C2H 4 M„ com- 
plex from the sum of the total energies of free molecular 
ethylene and of the M atom: E B (M) = [£; T (C 2 H 4 ) + 
nE T (M) - E T (C 2 ll4M n )}/n. According to the E B (M- 
atomic) values shown in the table, most TMs that we 
studied are able to bind relatively strongly to a C2H4 
molecule, except Cr and Zn. In table I, the variation of 
the TM binding energy with the number of the TM-3d 
electrons displays a behavior similar to what observed 
previously for the chemisorption of TMs on the surface 
of a single-walled carbon nanotube^i^. Namely, there 
exist two energy maxima between a minimum that oc- 
curs for the element with five ci-electrons. Table I also 
gives the binding energies with respect to bulk metal en- 
ergies (Eb (M-bulk)) while Table II also gives Eb with 



respect to metal dimer energies (Eb (M-dimer)). All 
E B (M-bulk) values are negative, indicating endothermic 
reactions. Apparently, metal atoms in vapor or some 
metal-precursors, instead of bulk metals, are preferred 
when synthesizing these complex structures. 

We next studied the H2 storage capacity of the metal- 
ethylene complex, by calculating the interaction between 
C2rl4M n and different number of H2 molecules. We con- 
sidered various configurations for the hydrogen absorp- 
tion on a metal center, as shown in Fig. 3. The first H2 
molecule absorbed may either be in molecular form (Fig. 
3(a)) or in dissociated form (Fig. 3(b)). For most transi- 
tion metals, it is possible to absorb more, up to five H2 
per M atom. Two of the many possible multiple H 2 ab- 
sorption configurations are shown in Fig. 3 (c) and (d). 
For Li, in both C2H4Li and C2H4Li2 complexes, each Li 
can bind to two H2, resulting in absorption capacity of 
10.3 wt % and 16.0 wt % respectively. The optimized 
configurations and structural parameters are shown in 
Fig. 4. 

The nature of the metal-H2 interaction is easy to un- 
derstand. For TMs, since the bonding orbitals are mainly 
between metal d- and hydrogen cr*-antibonding orbitals, 
the mechanism of this interesting interaction can be ex- 
plained by the Kubas interaction^. For Li, the metal- 
H 2 binding is mainly electrostatic. We summarize the 
average H2 binding energy for C2H4M in Table I. Note 
that the H2 binding energies differ slightly from those 
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FIG. 2: (Color online) Electronic density of states of C2H4, 
Li atom, and C2H4+Li complex. The isosurfaces of the molec- 
ular orbitals are also shown. The hybridization of the Li-2p 
state and the LUMO of C2H4 is apparent. See text for further 
explanation. 



(a) M+H 2 (c) MH 2 +3H 




(b) MH 2 (d) M+5H 2 




FIG. 3: (Color online) Various configurations that we con- 
sidered in this study, for the hydrogen absorption on a metal 
center of a C2H4M n complex: (a) one H2 absorbed molecu- 
larly; (b) H2 dissociating with two M-H bond formed; (c) two 
atomic H and three H2 molecules; (d) five H2 absorbed as 
molecules. Large and small balls represent M and H atoms, 
respectively. 



given in our earlier work for the C2H4M2 complexes with 
the sandwich structure. Nevertheless, in most cases, the 
H2 binding energies have the right order of magnitude 
for room temperature storage. Since the hydrogens are 
mainly absorbed molecularly, we also expect fast absorp- 
tion/desorption kinetics. 

In order to test their stability, we further studied the 
dynamic of the C2H4M n complexes by normal mode anal- 
ysis. We found no soft (i.e negative) mode, indicating 
that the complex structures are stable. Characteristic 
phonon modes are summarized in Table III, using Li 



(a) C 2 H 4 Li-2H 2 
(10.3 wt%) 



* * 

d(CC)=l,48 A 
d(LiC)= 1.99, 2.07 A 
d(HH)=0.78A 
d(LiH)=1.84, 1.92 A 
E B (H 2 ) per H 2 =0.24 eV 



(b) C 2 H 4 (Li-2H 2 ) 2 
(16.0 wt%) 



m * 

* * 

m 
* 

* 

1! 



FIG. 4: (Color online) Hydrogen absorption configurations 
on (a) C2H4Li and (b) C2H4Li2 complexes. Note that in both 
cases, each Li can bind two H2, resulting in high absorption 
capacities. Large, medium and small balls represent Li, C 
and H atoms, respectively. 



and Ti as examples. Our calculated mode frequencies 
for the C2H4 molecule agree very well with the experi- 
mental values^. Metal binding to C2H4 elongates and 
thus softens the C=C bond, resulting in lower stretch- 
ing mode frequencies. Also the softening of the CH2- 
torsion and CH2-bending modes is obvious. There are 
three main M-related vibrational modes. In two of these 
modes, M atoms vibrate parallel and perpendicular to 
the C=C bond. In the third mode, metal atoms vibrate 
perpendicular to the C2H4 plane. These three modes are 
unique for the C2H4M n complex and therefore should be 
present in any Raman/IR spectra of a successfully syn- 
thesized material. 

We also calculated the normal modes of C2H4M„ com- 
plexes absorbed with H2 and did not find any soft modes 
either, indicating that the configurations that we consid- 
ered indeed correspond to local energy minima. Among 
many vibrational modes, we note that the H2 stretch- 
ing mode is around 330-420 meV for the absorbed H 2 
molecules, significantly lower than ~540 meV for the 
free H2 molecule. Such a shift in the mode frequency 
would be the key feature that can be probed by Ra- 
man/IR measurement to confirm a successful synthesis 
of the structures predicted here. In the lower energy 
range, there are many M-H modes that are unique to the 
complexes. To manifest the M-H dynamics, we show in 
Fig. 5 the phonon density of states of O^^Ti^-Ha, com- 
plexes weighted by neutron cross-sections (note that H 
has much larger neutron scattering cross section than C 
and most metals) . These plots can provide a useful com- 
parison to experiments when trying to synthesize these 
materials. 
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TABLE III: Characteristic mode frequencies (meV) for C2H4, C2H4Li n and C2H4Ti n complexes. Experimental values for C2H4 
(from ref. 24) are also shown. Note that the metal-C2H4 binding significantly softens the C=C stretching, CHh-torsion and 
CHVbending modes. The three main M-modes give unique signatures for metal-C2H4 complexes. 



Mode/ Complex 


C2H4 C2H4, exp. 


C2H4Li 




C2H4Li2 


C 2 H 4 Ti 




C2H4Ti2 


C=C stretching 


202 201 


184 




172 


170 




167 


CHb-torsion 


128 127 


99 




49 


52 




56 


CH2-bending 


115-165 117-166 


84-145 




54-141 


94-140 




73-134 


M-vib, // C=C bond 




38 


40 (in 


-phase), 66 (out-of-phase) 


62 


15 (in 


-phase), 57 (out-of-phase) 


M-vib, _L C=C bond 




37 


22 (in 


-phase), 65 (out-of-phase) 


56 


22 (in 


-phase), 62 (out-of-phase) 


M-vib, 1 C2H4 plane 




40 


38 (in 


-phase), 74 (out-of-phase) 


63 


29 (in 


-phase), 48 (out-of-phase) 




11 ^LAL^hljL^^^^^ 



C,H,Ti+5rL 




800 
Energy (cm ) 



1600 



FIG. 5: (Color online) Simulated neutron inelastic spectrum 
for various C2H4Tin-Ha: configurations. Note that the M-H 
dynamics are unique and can be used as a probe to iden- 
tify these structures. Thus these plots can provide a useful 
comparison to experiments when trying to synthesize these 
materials. 
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IV. ACTIVATION ENERGIES AND REACTION 
PATHS BETWEEN DIFFERENT ISOMERS 

The C2H 4 M n and C2H 4 M rl -H. E complexes can have 
several isomeric structures. It is important to know the 
relative stabilities of these isomers and their implications 
for the hydrogen absorption properties. We thus studied 
the activation energies and reaction paths between dif- 
ferent isomers of C2H 4 M n complexes. Here we discuss 
representative results on M=Ti. 

We start with the C^ELiTi-l-H^ complex and con- 
sider two possible structural transitions, which lead to 
lower energy configurations through the dissociation of 
an H2 molecule over a Ti atom. In the first case, 
the H2 molecule dissociates on top of the Ti atom. 
C 2 H 4 (Ti-|-H2) and C 2 H 4 (TiH2) are the reactant and 
product, respectively. Their relaxed structures corre- 
spond to the first and last images shown in the top panel 



FIG. 6: (Color online) (a) The minimum energy path for the 
dissociation of the H2 molecule over the Ti atom complexed 
with C2H4. An energy barrier of ~ 0.25 eV is found for the 
dissociation. A total of 21 images were used in the NEB 
calculations, five of which are shown on the top. Marked 
circles in the potential plot are the points corresponding to 
these five images, (b) The activation energy plot for the for- 
mation of Titanol-molecule from C2H4Ti+H2 complex, indi- 
cating a very low barrier of «0.15 eV. Once the final product 
forms, the CCTi-bond angle is very soft, resulting in the zero- 
temperature structure shown in the inset. 



of Fig. 6(a). The calculated minimum energy path for 
this process gives ~0.25 eV barrier, which is small but 
still significant since the C2H4(Ti+H2) configuration cor- 
responds to a local energy minimum and possesses a H2 
binding energy of ~0.3 eV. In the second case, the H2 
molecule is first dissociated over Ti and then one of the 
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(a) 




(a) 




FIG. 7: (Color online) Hydrogen absorption on a titanol 
molecule, (a) One H2 binds molecularly to Ti, with a binding 
energy of 0.3 eV. (b) One H2 is dissociated, yielding TiH3 
structure. The corresponding binding energy is about 1.0 
eV/H2. (c) Five H2 bind as molecules to the titanol molecule 
with an average binding energy of 0.4 eV/H2. Large, medium 
and small balls represent Ti, C and H atoms, respectively. 



(a) 



(b) 





(c) 




0.2 0.4 0.6 0.8 

relative reaction path 

FIG. 8: (Color online) (a) Bare C2H4Ti2 dimer-perp com- 
plex, (b) The complex with seven H2 absorbed. Note that 
there also exists other stable configurations that not discussed 
here, (c) The minimum energy path for the transition from 
C2H4Ti2 sandwich configuration to dimer-perp configuration. 
The activation energy is about 0.55 eV. Note that regardless 
of which isomer of C2H4Ti2 we have, the resulting complex is 
able to bind multiple hydrogen molecules. 




0.2 0.4 0.6 Hi! 
relative reaction path 

FIG. 9: (Color online) The minimum energy paths for the 
transitions (a) from the C2H4Ti2 sandwich to the dissociated 
configuration and (b) from C2H4Ti2 dimer-perp to the disso- 
ciated C2H4 configuration, respectively. In both cases, there 
are large energy barriers on the order of 1.1 eV. 



H atoms goes to carbon, forming a CH3 group. The acti- 
vation energy plot for this process is shown in Fig. 6(b), 
indicating a very low barrier of only —0.15 eV. Once the 
product (i.e., the last image of the top panel of Fig. 6(b)) 
forms, the CCTi-bond angle is very soft, resulting in the 
zero-temperature structure shown in the inset, which has 
only 30 meV lower energy than the product. The final 
structure of the molecule (inset) is isostructural to the 
"ethanol" molecule and therefore we call it "titanol" . 

Since the titanol molecule is fairly easy to form, it is 
important to check if this new complex still possesses 
the high-capacity H2 absorption property. In Fig. 7, we 
show several stable hydrogen absorption configurations 
on a titanol molecule. With only one H2, it can be ab- 
sorbed molecularly (Fig. 7(a)) with a bind energy of 
0.3 eV or absorbed dissociatively (Fig. 7(b)), yielding 
a TiH3 structure, which has a binding energy of about 
1.0 eV/H2. We expect that the dissociation process may 
have a similar barrier to that found in Fig. 6 (a). Im- 
portantly, the titanol molecule can bind up to five H2 as 
molecules (Fig. 7(c)) with an average binding energy of 
-0.4 eV/H 2 . 

Next, we study the C2H4M2 dimer structures. For Ti, 
the dimer-perp structure (Fig. 8(a)) has lower total en- 
ergy than the isomeric sandwich structure (Fig. 1(b)) 
and dimer-par structure (Fig. 1(c)). Fig. 8(c) shows the 
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activation barrier for the transition from the sandwich 
configuration to the dimer-perp configuration. The acti- 
vation energy is about 0.55 eV. Shown in Fig. 8 (b) is one 
of the stable configurations that we identified for the hy- 
drogen absorption on the C2H4Ti2 dimer-perp structure. 
Apparently, regardless which isomer of C2H4Ti2 that we 
have, the complex is always able to bind multiple hydro- 
gen molecules. 

Finally, one may ask whether it is possible for the 
metal to catalyze and dissociate the C2H4 molecule (i.e., 
break the C=C bond), forming a more stable structure 
as shown in Fig. 1(f). Our calculations show that the 
activation energies for a sandwich to dissociated C2H4 
(Fig. 9(a)) and a dimer-perp to dissociated C2H4 con- 
figurations (Fig. 9(b)) are both large, ~1.1 eV. Thus it 
is very unlikely that the dissociation would happen un- 
der near ambient condition. Interestingly, we found that 
even the dissociated structure can still absorb multiple 
H2, in which case, the system is somewhat similar to a 
Ti metallocarbohedryne (met-car) cluste r 27 ! 28 . 



V. FINITE TEMPERATURE 
FIRST-PRINCIPLES MD SIMULATIONS 

In order to further test the stability of the C2H4M„-H :C 
complexes and the relative strength of different interac- 
tions (such as M-C2H4 versus M-H2 interactions) and 
to identify possible reaction paths, we have carried out 
extensive first-principles MD simulations in the micro- 
canonical ensemble (NVE) 22 i 23 . The system is first op- 
timized and then random initial velocities are generated 
to yield twice the target-temperature. When the system 
is in equilibrium, half of this energy goes to the poten- 
tial and therefore the final temperature oscillates around 
the target temperature. We note that due to the small 
atomic mass of some elements (e.g., Li and H) in our 
system, it is essential to use a small MD time step such 
as 0.5 fs. Furthermore, convergence criteria for energy 
at each MD iteration should be very accurate (we used 
10 -7 eV) in order to avoid total energy/temperature drift 
(i.e., change in the total energy/temperature as a func- 
tion of simulation time). Since we are studying an iso- 
lated molecular complex in free space, it is also important 
that we eliminate the six degrees of freedom (i.e. three 
rotations and three translations) of the molecule. When 
this is not done, we observed that the input tempera- 
ture goes to totally uniform translation/rotation of the 
molecules rather than populating the vibrational modes 
after 1-2 ps simulations. In our simulations, we fixed one 
of carbon atoms and then two components of position of 
the other carbon atom and one component of M atom 
position (which prevents the rotation of the molecule in 
the CCM-plane). In this way, the total degrees of free- 
dom allowed in our simulation are Np = 3 x (N — 6), 
as expected for an isolated molecule. The temperature 
of the system is defined as T(t) = J2i m-ivf /(2KbNp) 
where i runs over the atoms of the complex and fc^ is 
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FIG. 10: (Color online) First-principles MD results at 500K 
for the C2H4Li2+4H2 complex. Shown are the time evolu- 
tion of various quantities, including total energy (a) and tem- 
perature (b) of the system, C-C (c) and Li-C (d) bond dis- 
tances, and Li-C-C-Li torsion angle (e). The bottom panel (f) 
shows the distance between Li atom and the hydrogen cen- 
ter of mass, indicating successive desorption of H2 molecules 
along the simulation. 



Boltzman's constant. The relative fluctuation is of the 
order of 1/y/Np. We also note that since our system 
is very small (i.e. about a dozen of atoms), it is basi- 
cally a collection of a small number of harmonic oscil- 
lators and therefore temperature fluctuations are large. 
In fact trying to control system temperature through ve- 
locity scaling 2 ^^ at small time interval does not work 
and yields wrong results. The microcanonical ensemble 
is thus the best for our purpose and as we shall see below 
it works well provided that a small time step is used and 
the total energy/force calculations are accurate enough. 
Here we present representative results on M=Li, Ti as 
examples. 

Our MD results for C2H 4 Li 2 -l-4H2 at 500 K are sum- 
marized in Fig. 10. The constant of motion plot shows 
only 50 meV drift in total energy over 10 ps simulation 
time, which causes a small temperature drift. The C-C 
and Li-C distances, shown in Fig. 10(c) and (d) respec- 
tively, indicate that the bare C2H4Li2 molecule is stable 
at this temperature. The torsion angle Li-C-C-Li shows 
no sign of Li-dimer formation and oscillates around 180 ° . 
The bottom panel in Fig. 10 shows the distance between 
Li atoms and the center of mass of H2 molecules, indi- 
cating the successive release of hydrogen molecules from 
the system. The first H2 leaves the system around 400 
fs. The fluctuations in the distances become very large 
at 2000 fs, resulting from the release of another hydro- 
gen molecule. Around 8-10 ps, the other two hydrogen 
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molecules also leave the system. Even though with 10 ps 
MD simulations, it is not possible to get reliable temper- 
atures; the results are still very promising and suggest 
that the C2H4Li2 system can stay intact at 500 K while 
it releases four hydrogen molecules. 

We next studied the stability of C2H 4 Ti„ system. We 
performed MD simulations up to 10 ps on C2H4Ti2 
(sandwich), C 2 H 4 Ti+H 2 , C 2 H 4 (Ti+5H 2 )2 (sandwich), 
and the titanol molecule (CH 3 CH 2 TiH) at 300 K and 
500 K. In the simulations on the two metal sandwich 
systems, we did not observe any Ti-dimer formation. In 
the case for C2H4Ti+H2, we did observe the spill-over ef- 
fect, where the H2 is dissociated over Ti and C and then 
Ti moved away with one hydrogen atom attached to it. 
This is essentially the titanol formation process that we 
discussed in the previous section. 

The MD results for the C 2 H 4 (Ti+5H 2 )2 system at 500 
K are summarized in Fig. 11. During the 10 ps simu- 
lation time, both C-C and Ti-C bond distances oscillate 
around their equilibrium lengths without any indication 
of instability. Similarly, the Ti-C-C-Ti torsion angle also 
slowly oscillates around its equilibrium value of 180 ° and 
does not show any evidence for Ti-Ti dimer formation for 
which the torsion angle is supposed to be about 57 ° . Fig. 
11(d) shows the number of H2 molecules that are close 
to a Ti-atom (within a 2.2 A distance), showing that ini- 
tially two H2 molecules are released and then another H2 
molecule is released at around 2.4 ps. Above 6 ps, the 
number of H2 fluctuates indicating that the distances are 
going beyond 2.2 A more often. Probably if we had run 
the MD simulation further, we would lose the remaining 
H2 molecules that are attached to the Ti atoms. 

As a final example, in Fig. 12, we present results from 
a 10 ps MD run on titanol+5H2 molecules at 500 K. The 
C-C and Ti-C distances indicate that the bare titanol 
molecule is stable at this temperature. The C-C-Ti an- 
gle shown in Fig. 12(c) indicates that C-C-Ti bond angle 
is very soft, exhibiting large amplitude motion. Around 
5 ps, Ti actually goes to the middle of two carbon atoms, 
returning to our original C2H<4Ti like configuration. As 
we discussed in the previous section, these two configu- 
rations are almost degenerate. The last panel shows the 
number of H atoms that are within 2.2 A of the Ti atom. 
Three successive constructive desorptions of H2 molecule 
are evident. 

In summary, our MD results discussed above on dif- 
ferent systems indicate that the sandwich configuration 
of C2H4Ti2 is quite stable and can bind H2 molecules 
and then release them at elevated temperature. Simi- 
larly, C2H4Li2 MD results also suggest that Li is another 
promising option even though the strength of the inter- 
actions is at the low side. Finally, thanks to MD simula- 
tions, we discovered a new configuration, titanol, which 
is derived from the C2H4Ti+H2 system and capable of 
binding five H2 molecules and then releasing them at high 
temperature without breaking down its structure. 
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FIG. 11: (Color online) First-principles MD results at 500K 
for the C2H4(Ti+5H2)2 sandwich complex. Various quantities 
are shown, including C-C (a) and Ti-C (b) bond distances, 
and Ti-C-C-Ti torsion angle (c). The bottom panel (d) shows 
the number of H2 molecules that are within 2.2 A of Ti atoms. 
It indicates successive desorption of H2 molecules in the course 
of the simulation. 
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FIG. 12: (Color online) Various quantities (same as in Fig. 11) 
obtained from MD simulations of Titanol+5H2 system at 500 
K. 



VI. CONCLUSIONS 

Our conclusions are summarized as follows: 
1. We showed that the C=C bond in ethylene can 
mimic the double bond in other carbon structures like 
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Ceo, in terms of binding metal atoms and the hydro- 
gen absorption properties. The small system size of 
the M-ethylene complex allowed us to do very detailed 
studies such as long MD simulations and reaction path 
calculations, which were very difficult to perform other- 
wise. Most of the results that we found, such as H2- 
dissociation and titanol formation, should be valid for 
other Ti-decorated nanostructures. 

2. For light transition metals, we showed that the ini- 
tial H 2 absorption could be either molecular with binding 
energy of ~0.3 eV or it could be chemical by TiH 2 forma- 
tion with a binding energy of ~1.0 to 1.5 eV. However, 
there is a barrier of ~0.25 eV for this process. Since 
the molecular H2 has a binding energy of ~0.3 eV, the 
dissociation could not be observed. Indeed, in our MD 
simulations, we did not see conversion from Ti+H2 to 
TiH2. Instead, we discovered that there is a very low 
energy barrier for the simultaneous dissociation of H2 
and formation of CH bonding (similar to spillover effect) 
through the Ti atom. For the case of C2H4Ti+H2, this 
reaction yielded a new molecule which is isostructural to 
ethanol and can bind five hydrogen molecules with an 
average binding energy of ~0.4 eV. 

3. We showed that the sandwich configuration of 
C2H4M2 is quite stable for both transition metals and 
Li. There are high energy barriers for the transition to 
dimer-configurations. Our 10 ps MD simulations did not 
show any evidence for dimerization. 

4. From our results, it is clear that C2H4M„ system 
could have a very rich phase diagram with different con- 
figurations. However, for all the isomer configurations 
that we have investigated, the complex is always able to 



bind multiple hydrogen molecules with high absorption 
capacity. Hence, these results suggest that co-deposition 
of transition/lithium metals with small organic molecules 
into nanopores of low-density materials could be a very 
promising direction for discovering new materials with 
better storage properties. 

5. We note that there are many existing experimental 
studies of small organic molecules with transition metals 
in gas phase by mass spectroscopy. In these experiments, 
the metal atoms are obtained by laser evaporation of bulk 
metal and then condensed with mixture of Ar-and ethy- 
lene (or benzene) gas onto a cold substrate. In this way, 
it was possible to trap M X (C2H4) J/ types of complexes 
in an argon matrix and do spectroscopic experiments on 
them. We hope that our study will reenergize these stud- 
ies with the focus on hydrogen absorption properties of 
these systems. It may be possible to use H2 rather than 
Ar to prepare these clusters in an H2-matrix. Such stud- 
ies would be very important as a proof of concept and 
that should be the current emphasis. 
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